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Abstract. - Using the Quantum Inverse Scattering Method we construct 
an integrable Heisenberg-XXZ-model, or equivalently a model for spinless 
fermions with nearest-neighbour interaction, with defects. Each defect in- 
volves three sites with a fine tuning between nearest-neighbour and next- 
nearest-neighbour terms. We investigate the finite size corrections to the 
ground state energy and its dependence on an external flux as a function of 
a parameter u, characterizing the strength of the defects. For intermediate 
values of u, both quantities become very small, although the ground state 
wavefunction remains extended. 
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PACS: 71.27.+a — Strongly correlated electron systems. 
71.30.+h — Metal- insulator transitions. 
75.10.Jm — Quantized spin models. 



Introduction. - Three recent experiments have demonstrated that persistent currents, peri- 
odic in the magnetic flux, exist in mesoscopic metal PJ and semiconductor M rings at very 
low temperatures. Surprisingly, though the current is found to be small, of the order of 
~ evp/L for single rings (vp is the Fermi velocity, and L the circumference), it is still two 
orders of magnitude larger than expected theoretically, at least for the metal rings studied 
in Q. In the latter, the electron motion is diffusive, i.e. the elastic mean free path is much 
smaller than the circumference. While it is well established that the Coulomb interaction 
gives an important contribution to the current for a measurement on an ensemble of rings 
||, the interaction effect in single rings, is far from being understood theoretically. 
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In this article, we consider a one-dimensional, interacting model in the presence of a 
magnetic flux, or equivalently, with twisted boundary conditions. We introduce very special 
"defects" into the model describing spinless fermions with nearest-neighbour interaction. 
Despite this inhomogeneity, the model remains integrable and we present exact results for 
the finite size corrections to the ground state energy, and its dependence on the magnetic 
flux, as a function of a parameter v characterizing the strength of the defects. Clearly, our 
investigation does not provide an answer to the questions raised by the experiments (there, 
the number of transverse channels is much larger than one). Instead, our work is closely 
related to, and an extension of, various recent theoretical studies [4-9] of quantum coherence 
in strongly interacting electron systems. 

Construction of the model. - Using the Quantum Inverse Scattering Method (QISM), we 
construct our model from the 1Z and L matrices of the Heisenberg-XXZ-model on an inho- 
mogeneous lattice as, for example, described in The central equation of the QISM is the 
Yang-Baxter equation, which guarantees that a scattering process factorizes in two-particle 
scattering processes and does not depend on the order of these. In order to construct a 
model with defects, we allow that the local C n matrix depends, in addition to the spec- 
tral parameter A, on a parameter u n , C n (X) = £(A + u n ). The transfer matrix is given 
by T(A) = Tr]T^Li £(A + v n ), where M denotes the number of lattice sites. To include 
twisted boundary conditions, we multiply the Cm matrix of the Heisenberg-XXZ-ring with 
exp (i(fia z /2). 

The Hamiltonian is then given as the logarithmic derivative of the transfer matrix with 



respect to A, at a specific value fLOfl . In particular, in the special case in which all v n = 0, 
we obtain the usual XXZ-model, which can be transformed to a spinless fermion model 
by a Jordan- Wigner transformation. For a general set of parameters, {z/ n }, it is difficult 
to determine the Hamiltonian explicitly, with one exception, namely where there are no 
defects on neighbouring sites, i.e. u n u n+ i = for all n. This is the situation we study in 
the following. As an illustration, consider a vanishing nearest-neighbour interaction, and a 
single defect at the site n\ characterized by the parameter v. The resulting Hamiltonian is 
given by 

M 

H = H° + K^v) = - £ (c+c n+1 + c+ +1 c n ) + ?4(„) (1) 

n=l 

= (1 ~ cdh^) ( c »i-i c »i + c «i c «i+i) - tanh (^X-iC ni+ i + h.c. (2) 

where the {c+} and {c n } are the standard fermion creation and annihilation operators. The 
generalization to r defects is straightforward (assuming u n v n +i = 0), 

H = H° + iTH I ne (v ne ) (3) 
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where ne denotes the location of a defect with strength v ni . An illustration is given in Fig. [l|. 
The expression for the Hamiltonian in the presence of a finite nearest-neighbour interaction 
is more lengthy but similar in structure, i.e. a defect located at ri£ affects the lattice sites 
ri£ — 1, rig, and ri£ + 1 only [JQ]. 

Single defect, no interaction. - As is apparent from Eq. (fj), for v = 0, the Hamiltonian re- 
duces to 7i°, i.e. the standard single-band tight-binding model (the hopping matrix element 
is chosen to be unity). In the opposite limit, v = oo, the lattice site n± is cut out of the 
ring. As a result, the model represents free fermions on a ring of M — 1 sites, however, with 
an additional phase factor e lSl , 5i = ir/2, for the hopping between ri\ — 1 and ni + 1, plus 
one uncoupled site. We emphasize that the parameters, cosh _1 (z/) and tanh(z/), as well as 
the phase factor 5\ = n/2, are fine-tuned in the following sense: a generic impurity breaks 
translational invariance and lifts the degeneracies of the single-particle spectrum, which are 
found at = 0, ±7r. While our defects also break translational invariance, this symmetry 
is replaced by another, of not as clear physical origin. As a result we find that even when 
changing z/, no degeneracies are lifted — they only occur at different, //-dependent values of 
0. The corresponding //-dependent symmetry operators can be constructed [JTOj . 



The localization of electronic states is another, well established phenomenon, in one- 
dimensional disordered systems. In Fig. |[ we plot the squared modulus of the wave- 
function for the single-particle level lowest in energy. Clearly, for the integrable case, the 
wavefunction is extended though reduced at the defect. Allowing, however, the phase 5i to 
be different from tt/2, which corresponds to the non-integrable case, we find a drastically 
different behaviour with a clear localization of the wave-function near the defect. 

Several defects, finite interaction. - The results for a finite nearest-neighbour interaction, 
and in the presence of several defects, are obtained from the Bethe equations, which we derive 
from the algebraic Bethe ansatz, with the result 



cosh(Aj — irf) 
cosh(Aj + irj) 



M T tt cosh(Aj + u nt - if]) = ^ " sinh(Aj - \ k - 2irj) 
fJi cosh(Aj + v ni + irf) *Ji sinh(Aj - X k + 2ir]) ' 



Here, iV is the number of fermions and M and r denote the number of lattice sites and defects, 
respectively. The nearest-neighbour interaction, in units of the hopping matrix element, is 
parametrized as A = cos (2rj) (we consider — 1 < A < 1 only, A > corresponds to an 
attractive interaction). For the spin model and for the fermion model with an odd number 
of fermions, = (ir) represents periodic (antiperiodic) boundary conditions, and vice versa 
for an even number of fermions. The phase is directly related to the magnetic flux $, for 
odd iV according to the relation = 27r$/$ , $o = h/e. (For even N, = 27r$/$ - n.) 

The equations @ have a remarkable consequence: given a set of defect parameters {z/ n } 
the roots {A-,} of the Bethe equations do not depend on the distribution of the defects over 
the sites. This implies that the defects can be moved or permutated without any change 
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in the energy spectrum, which is quite different from the observation that in mesoscopic 
physics, thermodynamic and transport properties can vary considerably by moving only one 
impurity by a few lattice constants. 

In the following, we simplify our defect model further, by assuming that \v nt \ = u, 
and choosing an equal number of defects of strength +v and —v. With this choice, the 
ground state energy is an even function of the phase 0. We define x = lim^M^^r / M) 
to be the density of defects. Assuming in addition N = M/2, i.e. a half filled band, we 
calculate in a first step the ground state energy per lattice site in the thermodynamic limit, 

= limjvf^ooeM, &m = E M /M. Using the standard method HlQlJl2"ll, we find 



«. = - m r ,, A /°° (A) ,, A ,, d> - ( 5) 

j-00 cosh(A + ii]) cosh(A — ir\) 2 
p oo (A) = { 2(vr - 2 V ) cosh[7rA/(7r - 2r,)] (6) 
Poo (A) = (1 - x)p oo (X) + x [p^X + u) + p oo (A - u)]/2. (7) 

The last equation relates the density of roots for the infinite system, poo, to the corresponding 
density of roots of the homogeneous system, . For v = 0, clearly, p^ = p^, and we recover 
the results of |I2]] . On the other hand, for v — > 00, it follows from (|7|) that p^ = (1 — x)p oo ; 
as a direct consequence, the integral in Eq. (|5|) is simply multiplied by the factor (1 — x) 
compared to the homogeneous case, and e ro is easily calculated. Apparently, for all values 
of v, the ground state energy e ro is a linear function of the defect concentration (and it is 
independent of 0, i.e. the boundary conditions). 

The leading order finite size corrections can be determined with the help of the Wiener- 
Hopf technique, as described in [[|, with the result 

(tt 2 /6) sin (2^/(^-2^) / 30 2 ^ 1 
MW °°~ l-a; + xcosh[7r^/(7r-2r / )] { Anr] J M 2 ' [) 

which for x = is in agreement with |3[]. In contrast to e^, the finite size corrections 
depend in a nonlinear way on the defect density, and they are exponentially suppressed, 
~ exp [— ttv/(tt — 2t])], in the limit of large v. The finite size corrections, (e^ — e^M 2 , 
obtained by solving numerically the Bethe equations (^), are plotted in Fig. |3| as a function 
of the defect parameter v, for half filling (N = M/2), and A = 0.5, x = 0.2. The system size 
is varied from 10 2 to 10 4 . The numerical result is in perfect agreement with the analytical 
expression @, but we also observe that the limits M — > 00 and v — > 00 do not commute. 
Taking v — > 00 first, we obtain a system with r sites cut out; however, the occupancy of 
these sites still appears in the Hamiltonian, but as a good quantum number. The finite size 
corrections thus correspond to a system of M — r sites, i.e. they are enhanced in magnitude 
by the factor (1 — x)~ x = 5/4 as compared to a system with x = 0. For large values of u, the 
asymptotic result @ (applicable for v finite, M — > 00) is only reached for extremely large 
systems. 
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Finally, we consider the phase sensitivity of the ground state energy, similiar to Refs. [4- 
6], where the homogeneous system was studied. Here we choose M = 10 3 and x = 0.2, and 
we plot [ejif(7r) — e A/(0)]M 2 as a function of v in Fig. | for different fillings N/M. Close to 
half filling (N = 440, 490), the dependence on v is very similar to the result described above 
for the finite size corrections, i.e. the phase sensitivity is very small between v « 4 and v ~ 8. 
In this context, note that Eq. (||) implies that [eM^) — eA/(0)]/[eM(0) — e^] = — 371/(41]); 
but note also that different A's have been chosen in Figs. § and ^. 

The asymptotic (large v) results shown in Fig. || can be explained as follows: in this limit 
(and for the given parameters), 200 sites are cut out from the system of 1000 sites. But, as 
discussed above, the occupancy of these sites still enters into the Hamiltonian. The defect 
phases, 5i, are found to be ±2r/, and as a result, we find additional cusps in the energy-phase 
relation, for example, for rj slightly less than 7r/4, at <fi = ±4rj. Consequently, the number 
of discontinuities in the persistent current, I = —dEM/dQ, increases, implying a reduction 
of the phase average of the current over half a period, i.e. eui^) — cm(0), compared to its 
v = value. 

For low filling, N -C M, on the other hand, the interaction is unimportant and we may 
apply the free electron result, / = —(evp/L^/ir, which implies that E M {it) — E M (0) rs 
7r 2 N/M 2 for small N. (In our units, L = M and Hvp = 2 sin(ir N/M).) Comparing the 
limits v = and v — > oo, the phase sensitivity is thus approximately enhanced in the latter 
case by the factor (1 — x)~ 2 , which is apparent in Fig. |j. 

Conclusions. - We have presented exact results for the finite size corrections to the ground 
state energy and its flux sensitivity for an one-dimensional, interacting model in the pres- 
ence of defects. Through its construction, the model remains integrable. The defects are 
"magnetic" in the sense that a defect triangle, compare Fig. 0, encloses a finite magnetic 
flux. This means that time reversal symmetry is broken and, for general {z^ n }, the energy is 
not an even function of (p. Surprisingly, the energy spectrum is independent of the spatial 
distribution of defects, and we find neither level repulsion nor localized states, which are 
considered to be generic properties of a "real" impurity. We believe that the absence of 
these effects is strongly related to the integrability of the model. Nevertheless, the finite size 
corrections and the phase sensitivity, i.e. the persistent current, can become very small. 

We thank K.-H. Hock for interesting discussions and useful comments on the 
manuscript. P. Schmitteckert acknowledges financial support through a graduate student 
fellowship by the state of Bavaria. 
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List of Figures 

Figure [l|: Graphical representation of a defect at the site rig, compare Eq. (0). Note that 
for v — > (v — > oo) the dashed (dotted) lines representing the corresponding hopping 
contributions are effectively cut. 

Figure Squared modulus of the wavefunction of the lowest energy eigenstate for the 
non- interacting limit in the presence of a single defect (y — 1) at the site n\ = 50 
(M = 100). For the integrable case (Si = vr/2), the eigenstate is extended, though 
reduced at the defect, while a detuning of the defect phase (Si ^ 7r/2) leads to a 
localization of the wavefunction. 

Figure Finite size corrections to the ground state energy, (ejf — e^) M 2 , as a function 
of the defect parameter v for half filling, N = M/2. The system size is varied from 
M = 10 2 to M — 10 4 ; the results for large systems fit perfectly with the analytical 
expression, Eq. @ for = 0. The nearest-neighbour interaction is chosen to be 
A = 0.5, and the defect density is x = 0.2. 

Figure Phase sensitivity of the ground state energy, [6m(4> = tt) — cm(0)]M 2 , as a 
function of v for different fillings, i.e. numbers of fermions, iV (M = 10 3 , A = 0.2). 
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